Spin-spin coupling in electrostatically coupled quantum dots 
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We study the spin-spin coupling between two single-electron quantum dots due to the Coulomb 
and spin-orbit interactions, in the absence of tunneling between the dots. We find an anisotropic 
XY spin-spin interaction that is proportional to the Zeeman splitting produced by the external 
magnetic field. This interaction is studied both in the limit of weak and strong Coulomb repulsion 
with respect to the level spacing of the dot. The interaction is found to be a non-monotonic function 
of inter-dot distance do and external magnetic field, and, moreover, vanishes for some special values 
of do and/or magnetic field orientation. This mechanism thus provides a new way to generate 
and tune spin interaction between quantum dots. We propose a scheme to measure this spin-spin 
interaction based on the spin-relaxation-measurement technique. 
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I. INTRODUCTION 

The spin degree of freedom of electrons in semiconduc- 
tor nanostructures has attracted much interest recently. 
Besides its applications in spintronics 1 , which makes use 
of the spin in the same way as the charge in conven- 
tional devices, there has been considerable interest in the 
quantum coherence of spin in different materials and ge- 
ometries at the nanoscale 2 . The long coherence times 
measured some time ago in semiconductors 3 supported 
the idea of quantum computing with spins localized in 
quantum dots 4 ' 5 -a rapidly growing field now, where re- 
cent experimental progress has been remarkable 6-13 . 

Electron spins in semiconductor nanostructures, how- 
ever, are not decoupled from the charge degree of 
freedom, one of the primary reasons for this cou- 
pling being the spin-orbit interaction. This coupling 
leads to many interesting phenomena in the physics of 
semiconductors 1 ' 2 , an important one being the ability 
to control the spin of the electron with electric fields 
acting on its charge degrees of freedom 14-19 . In quan- 
tum dots the spin-orbit interaction manifests itself as a 
weak perturbation when compared with the confinement 
energy. The measure of smallness is given by the ra- 
tio between the dot radius A and the spin-orbit length, 
Xso — the distance over which an electron travels and 
thereby precesses by an angle n about the intrinsic 'mag- 
netic' field induced by the spin-orbit interaction. Despite 
its smallness the spin-orbit interaction is very impor- 
tant for the coherence of the spin dynamics. For exam- 
ple, spin-orbit interaction allows for coupling of the elec- 
tron spin to (bosonic) environments, such as phonons 20,21 
or particle- hole excitations in quantum point contacts 22 , 
which in turn causes relaxation and decoherence of the 
spins. Moreover, electrons being charged particles inter- 
act via the long range Coulomb forces with each other, 
even if they are confined to well-separated quantum dots 
with no overlap of their wave functions. Through this 
electrostatic coupling and in combination with the spin 
orbit interaction, the spins of two electrons located in 



different dots become coupled even in the absence of tun- 
neling between the dots. 

In this paper we provide a detailed analysis of such an 
effective spin-spin interaction for lateral quantum dots in 
a configuration as shown in Fig. 1. We will see that the 
origin of this interaction is the "tidal" effect each of the 
electrons produces on the charge distribution of the other 
electron via electrostatic forces. Because of the spin-orbit 
interaction, the electric dipole moment (as well as higher 
moments) in each dot couples to its electron spin. As 
a result, the two spins experience an interaction resem- 
bling a magnetic dipole-dipole interaction 23 with effec- 
tive magnetic moments which can be strongly enhanced 
by up to a factor of 10 3 compared to the Bohr magneton. 
The magnitude of the spin-spin coupling obtained via this 
spin-electric effect can be efficiently controlled and even 
completely suppressed by adjusting external parameters 
such as the magnetic field direction, strength, and inter- 
dot distance. 

The spin-spin interaction can, in principle, be used 
to perform two-qubit operations as required in the spin- 
based quantum computing scheme 4 , because it entangles 
spins and can easily be switched on and off. We note 
that a similar mechanism for spin interaction based on 
electrostatic coupling was studied very recently in Ref. 24 
for vertically coupled quantum dots, and in Ref. 25 for 
the special case of one-dimensional quantum dots formed 
in semiconducting nanowires. Besides the differences in 
geometry and dimensionality, both of these works treat 
only the case of weak Coulomb interaction (compared to 
the level spacing), while we treat here also the opposite 
limit of strong Coulomb interaction where new and inter- 
esting features emerge. In the limiting case of strongly 
elliptical dots we recover the one-dimensional results ob- 
tained in Ref. 25. 

We emphasize again that in the present study we 
exclude tunneling and thus the type of spin interac- 
tion studied in the following is fundamentally differ- 
ent from the Heisenberg exchange interaction for which 
the presence of electron tunneling between the dots is 
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crucial 26-29 . Similarly, the combined effect of Heisen- 
berg exchange interaction and spin-orbit coupling 30-34 
is also based on tunneling and should be carefully dis- 
tinguished from the spin-orbit effect studied here. We 
also note that the Heisenberg exchange coupling allows 
typically for much stronger spin-spin coupling than the 
electrostatically induced one. For instance, in GaAs dots 
the Heisenberg exchange can reach values on the order 
of 0.1 meV —1 meV, which, as we shall see, exceeds the 
electrostatically induced spin coupling by three to four 
orders of magnitude. Nevertheless, the electrostatic spin 
coupling can prove useful for cases where it is difficult 
to get sufficient wave function overlap (needed for large 
Heisenberg exchange), and, moreover, it is also impor- 
tant to understand the electrostatic spin-effects in detail 
in order to get control over possible interference effects 
between different types of spin coupling. This will be 
for example of importance for spin qubit applications in 
order to minimize spin decoherence and gate errors. 

Finally, in view of experimental tests we propose a 
scheme to measure the spin-spin interaction in a double 
dot setup with a nearby charge detector. We propose to 
combine the spin-measurement technique of Ref. 8 with 
the entangling property of the spin-spin interaction and 
present a gate pulsing sequence that enables one to access 
the coupling constant in the time domain by measuring 
the occupation probability of a Zeeman sublevel. 

The paper is organized as follows. In Sec. II we intro- 
duce our model for the quantum-dot system and its basic 
Hamiltonian. In Sec. Ill we derive the effective spin-spin 
coupling Hamiltonian within our model. In this section 
we analyze in detail the two limiting cases of weak and 
strong Coulomb interaction, and obtain various forms of 
the spin coupling as function of magnetic field strength 
and orientation, as well as of the interdot distance. The 
most general form of the spin coupling tensor is given 
in the Appendix. In Sec. IV, we propose a scheme to 
measure the spin-spin coupling. There, we also discuss 
the effect of the hyperfine interaction with the nuclear 
spins on the measurement scheme. Finally, in Sec. V we 
discuss our results and end with some conclusions. 



II. THE MODEL 

Our system consists of two electrons each of which is lo- 
calized in a quantum dot, and the two dots are separated 
from each other, without tunneling between them. The 
system is composed of two gate-defined quantum dots 
in a two-dimensional semiconductor layer (e.g. GaAs or 
In As) . A schematics of the system we consider is shown 
in Fig 1. 

We model the system by a harmonic confinement po- 
tential, which, for simplicity is assumed to be the same 
for both dots. Each dot is assumed to contain one 
electron with charge — e and spin S — (ft/2)<r, with 
c = "Vj cr z ) being the Pauli matrices. The model 




FIG. 1: The figure shows a sketch of the model system which 
consists of two identical quantum dots in the xy-plane, sepa- 
rated by distance ao (measured from dot-center to dot-center). 
Si denotes the spin of electron i = 1, 2, A is the dot radius, 
and B is the external magnetic field. The respective orbital 
wave functions of electron 1 and 2 are assumed to have no 
overlap (i.e. tunneling between the dots is excluded). The 
remaining purely electrostatic Coulomb interaction between 
the electron charges leads, via spin-orbit interaction, to an 
effective coupling between their spins. This spin-spin inter- 
action depends sensitively on the orientation of B, with no 
component along it, and is proportional to B 2 . 

Hamiltonian consisting of several terms reads 

H = Ho + H z + Hc + H S Oi (1) 

where Hq is the energy of the two electrons in the con- 
finement potentials 

Here, Pi = —ihd/dri + (e/c)A(ri) is the 2D kinetic mo- 
mentum of the i-th electron at position r,, m* the ef- 
fective mass, c the speed of light, U(n) = (m*/2)w 2 r 2 
the confinement potential for the z-th electron which is 
assumed to be harmonic and isotropic, and A(r,) is the 
electromagnetic vector potential. The strength of the 
confinement energy is given by the frequency c^o- The 
second term on the right-hand side of Eq. (1) is the Zee- 
man energy of the two electrons, 

Hz = ^g^BB ■ (<7! + <x 2 ). (3) 

The third term in Eq. (1) is the unscreened Coulomb 
interaction between the two electrons, 



K|ri-r 2 + oo|' 

where k is the dielectric constant of the material and an. 
is the geometric distance between the two dots, namely 
between the potential minima ('center') of the dots. With 
this choice, we measure the distance for each electron 
from its own dot center. The last term in Eq. (1) is 
the spin-orbit coupling which for strong z-confinement is 
given by 

H SO = \P(-p>i + pi°i) + a (p>i - Pv a i)i ( 5 ) 

i=l,2 
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being the sum of the Dresselhaus term 35 {j3) coming 
from bulk inversion asymmetry and the Rashba term 36 
(a) coming from structure inversion asymmetry. We as- 
sume the same coefficients (3 and a for both dots. It 
is convenient to work with center-of-mass and relative 
coordinates 37 , as the Coulomb interaction couples only 
to the relative ones and the solution of the center-of- 
mass part is straightforward 38 ' 39 . This then involves 
the standard substitutions M — 2m*, m = m*/2, 
R = (n + r 2 )/2, r = n - r 2 , and P = p x + p 2 and 
P = (Pi -P2)/2. 

III. SPIN-SPIN COUPLING 

We now turn our attention to the spin-orbit interac- 
tion. As was shown in Ref. 21, the spin-orbit coupling 
gives non-zero first order effects only if a magnetic field is 
present, as a consequence of the Kramers degeneracy. In 
order to describe the effective first order spin-orbit term 
in the presence of a magnetic field we make use of the 
Schricffer- Wolff (unitary) transformation 19,21 

e s (H d + H so )e- s =H d + AH, (6) 

where S = — is chosen such that PAH = AH, with 
the projector operator P satisfying PA = J2 n ^™l n )( n l 
VA, and H d \n) = E n \n). The Hamiltonian H d = H +H c 
(or H d = H R + H r + H c in center-of-mass and relative 
coordinates). The Hamiltonian H = H d + AH is diago- 
nal in the basis of H d and has the same energy spectrum 
as the Hamiltonian H = H d + Hso- In first order of the 
spin-orbit interaction Hso the transformation generator 
becomes S = (1 — P)^ 1 Hso, where L d is the dot Liou- 
villcan, L d A = [H R + H r + H C ,A], VA Evaluating this 
expression explicitly we obtain 

S=(l-P) (7) 

i=l,2 

with^i t 2 = (t/i,2/A_|_, Zi^/A-, 0). In this lowest order the 
transformed Hamiltonian H becomes 

H = H R + H r + H z + H c + H§ , (8) 

with H§ Q = H§^ + H§£ where 

H z s x 6 2 = [S,^ 2 ] = Ez(l-P)[lx{t R ±t r /2)]-v li2 . (9) 

In Eq. (9) I — B/B is the magnetic field direction vec- 
tor, Ez = 9^bB is the Zeeman energy and the vec- 
tors £ R and £ r are given by £ R = (Y'/A_,X'/\ + ,0) 
and £ r = (y'/A_, x'/X + , 0), respectively. The new co- 
ordinates correspond to a rotation by an angle 7r/4 — 7 
with respect to the coordinate frame in which the direc- 
tion of the ao-vector is associated with the x-axis in the 
XY(xy) plane so that the final expressions have the sim- 
plest form 21 . Here, 7 is the angle between the xy frame in 
Fig. 1 and the normal axes of the crystal. The spin-orbit 



lengths X± are given in the form 1/A± = m*{[3 ± a). 
We are now in a position to derive the coupling be- 
tween the two spins. This is achieved by performing a 
second Schricffer- Wolff (SW2) transformation which ex- 
cludes the first order contribution in spin-orbit interac- 
tion with no diagonal matrix elements. The new trans- 
formed Hamiltonian has the form H e g = e T He~ T , with 
T given by 

T=(l-P)(L d + L z y 1 Hl . (10) 

We assume now that the Zeeman energy is smaller than 
the orbital confining energy, Ez <C huj , which is usu- 
ally the case for electrons in quantum dots, such that we 
can neglect in Eq. (10) the Zeeman Liouvillean Lz (for 
spin-orbit effects due to level crossing see Ref. 40). In 
second-order in spin-orbit coupling the effective Hamil- 
tonian H e ff becomes 

H e ff = H d + H z + -[L^HgQ^go]. (11) 

The last term in Eq. (11) contains the desired spin-spin 
coupling between the two spins. However, besides this 
interaction it also contains some self-interaction terms 
which rcnormalizcs only the Zeeman splitting. We will 
not study those terms since they are of no practical in- 
terest in the case of identical dots. We consider a general 
magnetic field B — (cos $ sin 8, sin $ sin 9, cos 9), where 
9 is the angle between the magnetic field and the z-axis 
perpendicular to the 2DEG plane and $ the angle be- 
tween the in-plane component of the magnetic field and 
the x-direction (Fig. 1). The interaction between the 
two spins has the most general form 

Hs = \Y} L d lH lb,H Z so ], i,j = l,2. (12) 

The spin Hamiltonian is obtained by averaging over the 
orbital ground state, H s = (0|i? s |0). We then obtain 

H s =<Tl-I(T2, (13) 

where 

W ab = E 2 z (0\[(l x V£i)„, (/ x 6)b]|0), a,b = x,y, z. 

(14) 

We note that there is no component of the spin along the 
magnetic field direction as a consequence of the vector 
product in the tensor M. By diagonalizing the above 
tensor, we obtain for the Hamiltonian H s the reduced 
expression 

- Jx<t14 + Jy°l a l ( 15 ) 

where the couplings Jj^ depend on the magnetic field ori- 
entation and on the functions C ai b 2 = (0| [L d ai, b 2 ] |0), 
with a, b — x,y (for explicit expressions see the Ap- 
pendix). Thus, the effective spin-spin interaction is 
highly anisotropic, and, in general, of the AF-type. We 
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note in particular that for an in-plane magnetic field 
(9 = 7r/2), the spin Hamiltonian reduces to the Ising 
Hamiltonian, H s = JyO^a 2 (in a transverse magnetic 
held). Next, we rewrite H s in terms of raising/lowering 
spin operators a± 



where the operators Spi, i—1,2, are the charge density 
distribution modifications in each dot as a consequence 
of the spin-orbit interaction. They are defined as 



Sp l =p i -p°, i = 1,2, 



(20) 



H s = Jeg{o\o 2 _ + o\a\) + J' eff (a 1 _a 2 _ + a^crj.), (16) with Pi being the charge density operator in the absence 



with J eff = (1/2)(J S + Jy) and J' eff = (1/2)(J 5 - J 5 ). 
We recall now that the full spin Hamiltonian includes 
the Zeeman energy, given in Eq. (3), which leads to a 
large energy gap with 2Ez S> Jz,y- We will find below 
that typically 



Jx,y Ez ( A 
Ez hiv \Aso 



< 1 



(17) 



under our assumption that Ez <C fruJo and A <C Xso- 
As a consequence, we can neglect in Eq. (16) the terms 
proportional to J' e g since they cause transitions between 
different Zeeman levels of the total spin. The relevant 
spin-spin interaction, , which acts only within the 
S — T subspace, becomes then 



Hf = J eff (°W- 



(18) 



Thus, we are left with the task of calculating the coupling 
strengths and J e g. Because of the Coulomb term, 
Eq. (4), this cannot be done exactly and some approx- 
imations are required. They will depend on the ratio 8 
between the Coulomb interaction strength, e 2 /na , and 
the orbital level spacing, fiw , giving 8 — (e 2 / nao) / fhoo = 
(A/as) • (A/ao), with A = yJh/m*u>o being the dot radius 
and cib = ti 2 n/m*e 2 - the Bohr radius in the material. 
In other words, the parameter 8 will dictate the physics 
of the system, and from now on we will speak of the ratio 
\/a B as being the Coulomb interaction strength (repre- 
senting in fact the 'true' Coulomb strength for touching 
dots). For making the following analysis more transpar- 
ent we focus on the case with only Rashba spin-orbit 
coupling (A_ = A + = Ago)- The generalization to the 
case with both Rashba and Dresselhaus terms present 
is straightforward, but at the cost of more complicated 
expressions (see Appendix). 



A. Weak Coulomb coupling - 

One interesting case is met when { < 1, such that 
the Coulomb interaction can be treated as a perturba- 
tion compared to the orbital level spacing. In this case, 
one can retain only the first order contribution from the 
Coulomb interaction, which translates into the approxi- 
mation L^ 1 L^ 1 — Lq 1 LcLq 1 . Making use of this and 
after some algebra we obtain for the spin-spin coupling 
the following expression 



H, 



-I 



dr\dr2- 



8p\8p 2 



K\ri -r 2 + o 



(19) 



of spin-orbit interaction and pi = e T "p®e~ T ° the one in 
the presence of spin-orbit interaction, with To = L^H^q 
for the present approximation. From Eq. (19) we see 
that the spin interaction results from a Coulomb-type of 
coupling between two charge density distributions which 
themselves depend on spin. 

Let us now analyze in more detail Eq. (19). The first 
task is to find 8pt, for i = 1, 2, namely the spin-orbit in- 
duced charge distribution or the spin-dependent charge 
distributions for each dot. In order to do this, we give 
first some important relations valid in the case of har- 
monic confining potential, relations which are used in 
the following for the derivation of the main results 



L7 1 



eB z 

hni*w 2 V K + ~ y ' 



o Vi = ~ 



hm*uj 2 



Pi- 



- X i 



L o 1 Pi = — 'fi- 



lm 



(21) 



(22) 



(23) 



Making use of the relations Eqs. (21-23) and within the 
first order of spin-orbit coupling, i.e. Spi w [T ,p®], we 
obtain 



2E z e 



m*X 2 uj 2 Xso 



pot [ cos 6(yi cos $ + Xi sin <&)<t* 
{yiSm$-XiCOS$)ai], (24) 



with p° being the bare charge density in the dot cor- 
responding to the ground state and which assumes the 
well-known form for harmonic potentials 



P?(r) = 



(g? + Vi) 
A 2 



(25) 



We note that when there exist a perpendicular com- 
ponent of the magnetic field, the dot radius is renor- 
malized d ue to t he orbital effect of the magnetic field 
A — > A/Vl + r 2 , with r = uj /2lo c (uj c — eB z /m*c 7 
B z = Bcos6). However, we will still refer to A as be- 
ing the dot radius, with the appropriate expression de- 
pending on the magnetic field orientation. We could now 
insert the expression Eq. (24) for Spi in Eq. (19) and 
compute directly the spin Hamiltonian. However, work- 
ing with the Coulomb potential, it is more convenient 
to work with the center-of-mass and relative coordinates 
and for simplicity the x-axis along the inter-dot direction 
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ao. Assuming for simplicity a perpendicular magnetic 
field, the spin Hamiltonian H s takes the form 



H, 



4£ 2 e 2 



JJ dr dR Po (r)po(R) 
(X 2 - X 2 /A)alal + {Y 2 -y 2 /A)ay y 



m* 2 X^\% 



K Vy 2 + + a o) 2 



(26) 



with the electronic densities Po{r) = 
(2/7rA 2 )cxp(-r 2 /2A 2 ) and p (R) 

(l/27rA 2 )exp(-2i? 2 /A 2 ). In Eq. (26) there are no 
mixed terms like a].a 2 since those terms vanish because 
of the odd symmetry of the integrands in the case of 
harmonic confinement, which reflects inversion symme- 
try. The limit of in-plane magnetic field is obtained very 
easy from Eq. (26) by substituting the denominator 
with [{X 2 - x 2 /4) cos 2 $ + (Y 2 - x 2 /4) sin 2 §\o\o 2 y. [For 
general field orientation the expression for H s is more 
complicated (see Appendix).] In order to make the 
following analysis more transparent, we introduce the 
dimcnsionlcss coordinates r — > r/X and R — > R/X. The 
integration over the center-of-mass coordinates is now 
straightforward and the reduced expression for the spin 
Hamiltonian becomes 



H. 



El 



{AEZalo* + A££a> 2 ) , (27) 



for a perpendicular magnetic field and 



H„ = 



E 2 



m* 2 X^X 2 so 



(AE* sin 2 $ + AE y c cos 2 $) a\o\, 



(28) 

for an in-plane magnetic field oriented at an angle $ with 
respect to the inter-dot distance vector a . The energy 
differences AE^ V are given by 



AE* c = ^f dr Po (r) _L 

«A 2 J Jyi + ( x 



AEl 



— I 

nX 2 J 



drp (r) 



yV + (x + a /A) 2 ' 



vV + (x + a /A) 2 ' 



(29) 



(30) 



The ground state and the first excited states of the dots 
in relative coordinates give rise to different charge dis- 
tributions (po ; Pixi and Piy, respectively), and thus to 
different potential energies seen by a test charge at a dis- 
tance ao (along x) away from the center of the charge 
distribution (in relative coordinates). AEq V are the dif- 
ferences between these potential energies. 

Before studying the distance dependence of the spin 
Hamiltonian H s (determined by AEq V ) in the entire 
range of distances, it is instructive to see how the expres- 
sion Eq. (19) behaves in the large distance limit, ao 3> A, 
and to make some comparison with the magnetic dipolar 
interaction in vacuum 2 . We perform a multipole expan- 
sion of the Hamiltonian in Eq. (19). The first non-zero 



contribution takes the form of a dipole-dipole interaction 
between two spin-dependent electric dipoles, or phrased 
differently, the interaction between two charge-induced 
magnetic dipoles 

mi • m 2 - 3(mi • n a )(m 2 ■ n a ) 
Hs M ^ (31) 

with the dipolc moments rrii given by 

m l = Tr orb [8p l r i ] = pcr l ,i = 1,2. (32) 

Here, the trace is taken over the orbital degrees of free- 
dom with n a = ao /ao and p being the tensor correspond- 
ing to an effective spin-orbit induced magneton 



eE z 



cos9 (T 
cosd | . (33) 



m*uj 2 X S o { o sin g o 

We see from Eq. (33) that the tensor p, depends on 
the magnetic field orientation with respect to the 2DEG 
and that it is also anisotropic, in contrast to the usual 
isotropic Bohr magneton ps = eh/2m e c (m e is the mass 
of the free electron and c the speed of light). We note 
that the z-component of the induced magnetic moment 
(with z\\B) vanishes, i.e. m = (m x ,m y ,0). Let us quan- 
tify the strength of p, by the norm p| | = (1/3) M^, 



i.e. 



eE z 



3 m*u^X S o 



a/1 + cos 2 6. 



(34) 



We compare now \\p\\ with ps- First of all, we note 
that | \p\ | vanishes when there is no Zeeman splitting. 
However, for finite magnetic fields, \\p,\\ can exceed ps 
by many orders of magnitude in the case of quantum 
dots. To give an estimate, we assume TvjUq w 0.5 meV, 
E z « 0.05 mcV (B f» 2T) and to* = 0.067m e , A so ~ 
10~ 6 m for GaAs quantum dots. With these values, and 
taking 6 = (perpendicular magnetic field) we obtain 



4 E z m e 



p B 3 fraj m* ujqXso 



10 d . 



(35) 



We describe now in more detail the limit of large dis- 
tance between the dots. From Eqs. (27) and (26) we find 
for a ^> X 

H s = J(<ry y - 2al* 2 x ), (36) 

for a perpendicular magnetic field, and 

H s = J(cos 2 $ - 2 sin 2 <5>)o\(jI (37) 

for an in-plane magnetic field, with the coupling strength 
J having the form 



J 



o"so a o 



(38) 
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From Eq. (38) we see a large distance decay ~ 3 , i.e. 
a long range type behavior. We note also that the large 
distance result in Eq. (38) does not depend anymore on 
the orbital effect of the magnetic field. Working instead 
with the effective Hamiltonian defined in Eq. (f 8), the ef- 
fective coupling strength J e g for arbitrary magnetic field 
is given by 

Jeff = - ^ [l + sm2 0(1 - 3 sin2 $ )] ■ ( 39 ) 

We note that the spin coupling can range from maximally 
ferromagnetic interaction with J e g = — J at 9 = n/2 and 
<f> = to maximally antifcrromagnetic interaction with 
J e ff = J/2 at 9 = 7r/2 and <E> = 7r/2, and with J e g 
passing through zero for certain angles. 

Next, we consider the case of arbitrary distance but 
still with small Coulomb interaction strength as defined 
at the beginning of this section. Considering defined 
in Eq. (18) we obtain 

J eff = E Z — |f G(ao/A,M), (40) 

where 

G(a /A, 6, $) = ^ [(cos 2 9 cos 2 $ + sin 2 $)AE y c 

+ (cos 2 9 sin 2 $ + cos 2 $) A£§] . (41) 

The function G(a /A, 0, $) is plotted in Fig. 2 for dif- 
ferent angles 9 and <&. As for the large distance limit in 
Eq. (39), a similar but more complicated ferromagnetic- 
antiferromagnetic crossover behavior occurs as a function 
of the field orientation. However, in this case this behav- 
ior can also be induced by changing the distance between 
the dots a,Q (see Fig. 2). 

Eq. (40) suggests that the condition S <C 1 
is too restrictive. Instead, the weaker condition 
(A/as)(A/ao) 3 <C 1 is sufficient for the approximation 
to be valid. Fig. 2 shows a breakdown of the dipolar 
approximation (i.e. of the <z^~ 3 decay), occurring at a dot 
separation ao / A ~ 2 for perpendicular magnetic fields 
(9 = 0), and also a cancellation of this interaction for 
some given distance, which is around a /X w 1.8. This 
shows that the sum of the two electrostatic energy dif- 
ferences AEq + AEg has a non-monotonic behavior as 
a function of the distance ao- Actually, only AEq is 
non-monotonic, whereas AEq has a positive value which 
decreases with ao, as can be seen from Fig. 2. If an 
in-plane magnetic field is applied along y ($ = 0) or x 
(<I> = 7r/2) direction, we obtain a dependence only ei- 
ther on AEq or on AEq. Accordingly, G(ao/X) will 
be larger in some parameter range as compared with 
the case of perpendicular fields, see Fig. 2. At this 
point it is instructive to consider numerical estimates for 
the coupling strength J e g. For this we consider GaAs 
quantum dots for which we assume \so ~ 10~ 6 m, 
to* = 0.067 TO e , g = —0.44, k — 13, and also con- 
sider huj w 0.5 meV, E z = 0.05 meV (B « 2 T) and 
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FIG. 2: The function G occurring in Eq. (41) plotted as a 
function of the geometric distance ao between the dot centers 
scaled by the dot radius A for different magnetic field orienta- 
tions. The dashed line represents the dipolar approximation 
of G for a perpendicular magnetic field (9 — 0) which scales 
like a^ 3 . 



ao = 5 • 10~ 7 m. These estimations lead to a coupling 
strength J w 10~ 10 eV, which lead to a time dynamics 
of the order of 10~ 5 s. If this time scale is longer than 
the decoherence times in GaAs quantum dots, the system 
will be insensitive to the coherent dynamics induced by 
the coupling J e ff- Shorter time scales are obtained for 
materials with larger spin-orbit coupling such as In As. 
The spin-orbit length X$o m this material is comparable 
with a typical dot size of about 100 nm. Even though 
our perturbative approach starts to get unreliable in this 
case, it still can provide a rough estimate for the coupling 
strength. For InAs we have to* = 0.023to c , g = 14.8, 
k = 13, Xso ~ lOOnm, and we choose Hlvq f=a lmeV, 
Ez ~ O.lmeV and a /A w 3. With those values we ob- 
tain for the coupling J e g k, 10~ 7 eV, which corresponds 
to a switching time of about ~ 50ns for a swap of the 
spin states of electron one and electron two. This time 
scale for the spin dynamics is shorter than the expected 
spin decoherence time in such quantum dots. Thus, this 
interaction mechanism provides a useful way for the dy- 
namical control of the spin-spin coupling. As discussed 
before, for an in-plane magnetic field the coupling con- 
stant could even be higher, depending on the angle of the 
magnetic field with respect to the inter-dot axis. 



B. Elliptical dots with 5 < 1 

We briefly generalize the previous results to elliptical 
dot shapes. This will also allow us to study the one- 
dimensional limit and recover previous results obtained 
for one-dimensional nanowires 25 . We consider elliptical 
dots which are characterized by the frequencies ojq x and 
tooy corresponding to the x and y directions, respectively. 
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In this case, Eq. (38) is replaced by 

e 2 E 2 ( AE£ 2 2 x 2 \ ,. * 

where the electrostatic energies AEg' y become now 

"li Vy + O + ao/Ai) 2 



an Ising interaction which, together with single qubit ro- 
tations, can be used 4 to efficiently perform CNOT gate 
operations between two qubits. We finally note that in 
one dimensions the Rashba interaction can be treated 
exactly, leading to a renormalization of the (/-factor 25 , 
g — > .gcxp (— A 2 /A| ). This exact treatment is no longer 
possible in the 2D case considered here, except for the 
special case 42 when a = ±f3 and the problem becomes 
effectively ID 42 . 



f 



i-y 2 



A I:'' = — / dr Po (r)— 

° «W Vy 2 + (^ + ao/A 2 ) 2 



(44) 



with the charge density distribution function, expressed 
in relative coordinates, 

Po{ r) = ^-e-* 2 /^-y 2 /^l (45) 

For perpendicular magnetic fields and elliptical dots the 
lengths A 1,2 are given by 41 



Ai, : 



/4fi(n + l) 



(46) 



where n = m\m<2.W\WiU 2 l B 2 , Wi,2 = VA ± B/2, and 
toi, 2 = 2_B/(C±cj 2 + _B) with the explicit expressions for 
A, B, and C 



A 2 2,2 



(47) 



B 



c 



J 0y 



J 0x 



u 0y 



^l v , (48) 



(49) 



Taking now also the limit of strongly elliptical dots, i.e. 
ojQy » ojo x ,uj c , we see that this is equivalent to keeping 
only one component of the spin interaction, namely the 
<J x a 2 part, and that the orbital effect of the magnetic 
field drops out. The resulting Hamiltonian then becomes 



H„ = 



e 2 E 2 AE% 



c 



J 0x^ ^SO 



ial + Oii^/LJoy) 4 ). (50) 



Considering now the large distance limit, a 3> A, anal- 
ogously to Eq. (36), our result reduces formally to the 
one in Ref. 25, i.e. 



H„ = -2 



E *' 2 alal + Oii^/cJoyf)- (51) 



Km* 2 w$ x \ 2 so al 



The above expression, Eq. (51), could be obtained di- 
rectly from Eq. (57) since within the considered limit 
there is no orbital effect of the perpendicular magnetic 
field on the spin-spin interaction. We note that in this 
limit the resulting spin-spin coupling takes the form of 



C. Strong Coulomb coupling - S > 1 

We turn now to the more involved case of strong 
Coulomb interaction strength, 5 > 1, which cannot be 
treated perturbatively. However, some approximations 
are still possible and we will explore two of them in the 
following section. The first approximation consists in re- 
ducing the two electron system to two classical point- 
charge particles. The classical equilibrium condition will 
be obtained by minimizing the total potential energy of 
the two particles. By doing this, the motion of the elec- 
trons will take place around the new equilibrium posi- 
tions obtained from the equation 



a 2 (a — ao) = 2X 4 /as, 



(52) 



where oo is the initial geometric distance and a the ef- 
fective distance between the electrons in classical equilib- 
rium. However, we are interested in the motion around 
the equilibrium position, which means that for small de- 
viations, we may substitute the full Coulomb interaction 
with an effective one, remembering that r = n — r 2 , 



K\r + a \ 



2 k a- 



(3(n„ • r) 2 - r 2 ). 



(53) 



We note that the coordinates are measured now with 
respect to the new equilibrium position. Within this 
approximation, the relative Hamiltonian H r is replaced 
with the new, renormalized one, H r 



i 



i 



(54) 



with the definition u x , y — b x ^ y u; and b x>y given by the 



expressions 



b x = y/l + 4(\/a B )(\/a)3, b y = v /l-2(A/a B )(A/a)3. 

(55) 

We see that this approximation leads only to a renormal- 
ization of the orbital frequencies of the relative Hamilto- 
nian. We proceed now to develop a second alternative for 
treating the Coulomb interaction, namely a variational 
method based on the same picture of classical equilib- 
rium. This consists in substituting the full Coulomb in- 
teraction term with the same type of expression like in 
Eq. (52), but with the effective distance obtained from 
the variational ansatz. For this we minimize the expec- 
tation value of the orbital Hamiltonian, H r + He, in the 
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FIG. 3: a) The effective distance a as a function of the 
Coulomb interaction strength A/as for ao/A = 2. The full 
line represents the variational result from Eq. (56). The 
dashed line corresponds to the one obtained from the classi- 
cal equilibrium solution of Eq. (52). b) Effective distance a/ A 
as a function of the geometrical one ao/A for A/as = 3. All 
distances are scaled with the dot radius A. The dotted line is 
ao/A which is shown for comparison. 
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FIG. 4: The logarithm of infidelity 1 - F from Eq. (61) as 
a function of the geometric distance ao/A (scaled with the 
dot radius A) for two different values of Coulomb interaction 
strength A/as- 



ground state of the effective relative Hamiltonian H r [Eq. 
(54)] with respect to the effective distance a. This leads 
to the following equation 



— (ipo\H r + H c \ip ) =0, 



(56) 



where \ipo) is the ground state belonging to H ri i.e. 
H r \i>o) = -EolV'o): w ith E the ground state energy. Since 
we are dealing with harmonic oscillators, those wave 
functions are known. However, Eq. (56) for the effective 
distance a can be solved only numerically. We plot in 
Fig. 3 the results obtained for the effective distance a as 
a function of different parameters in both cases, namely 
the variational result from Eq. (56) and also the result 
obtained from the classical equilibrium condition in Eq. 
(52). We see in Fig. 3 that there is very good agree- 
ment between the two approaches in a wide parameter 
range and moreover, that the effective distance within 
the variational approach is larger then the one obtained 
from classical equilibrium, that means a lower ground 
state energy. We note that a perpendicular magnetic 



field practically does not change the curves in Fig. 3 
(not shown) on a large range of magnetic field strengths 
(0 < oj c < 3luq), which means that the effective distance 
is to a very good approximation independent of the ap- 
plied magnetic field. In order to verify the accuracy of 
our variational method, we checked also the numerical fi- 
delity, defined as the overlap of the wave functions in the 
variational case with the exact (almost, in the sense of 
perturbation theory) wave function. Although the prob- 
lem contains no small parameter, we can still define small 
matrix elements compared with level spacing in a numer- 
ical sense. For this we write the full relative Hamiltonian 
in the following way 



(57) 



with the effective Coulomb interaction V (see Eq. (53)) 
expressed now also in terms of the new equilibrium coor- 
dinates (introduced after Eq. (53)), 



V = 



Kir + a\ 



2kci' 



(3(n a0 • r) 2 - r 2 ) + 



2 

e x 



2 ' 



(58) 



We show now that this term leads to small matrix ele- 
ments so that indeed \Vo n \ "C \E n — Eq\, where the en- 
ergies E n and E are the n-eigenvalue and ground-state 
energy of the Hamiltonian H r , respectively. To see this 
numerically we introduce the fidelity 



F 



(^c#o) 



(59) 



where \ipo) and \ipo) are the ground state wave functions 
of the full Hamiltonian H r and H r , respectively. We now 
estimate the fidelity F by using perturbation theory to 
find the true ground-state wave function |^o) from the 
effective one l^o) 



\ipo) = \ipo) 



—\tpn) + -.-, (60) 



n=l 



E„. — En 



where we retain only terms to first order in V. Taking 
into account Eq. (60) we obtain the infidelity, 1 — F, 
namely the deviation of the true ground state wave func- 
tion from the effective one 



E 

n=l 



E„. — En 



(61) 



We plot in Fig. 4 the infidelity 1 — F as a function of 
the effective distance a for fixed Coulomb strength, X/as- 
We see that the infidelity takes very small values (1 — F < 
10~ 2 ) on the considered range, for two different Coulomb 
strengths, which shows that our variational approach is 
very accurate. 

We can now evaluate the spin-spin interaction within 
this approximation. Since we are now dealing with har- 
monic potentials only, the problem of finding J^.y from 
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FIG. 5: The function / from Eq. (67) as a function of the 
dimensionless geometric distance do/A. a) The case of per- 
pendicular magnetic (9 — 0) field for three different Coulomb 
strength parameters A/as. b) The case of in-plane magnetic 
field (6 = 7r/2), for two values of the Coulomb strength A/as, 
for three angles between the inter-dot distance vector ao and 
the magnetic field. The groups of lines 1, 2, and 3 correspond 
to $ = 0, $ = 7r/4, and <3> = 7r/2, respectively. 



Eq. (15) becomes straightforward. However, in the 
derivation of the spin Hamiltonian we need again some 
relations, similar to Eq. (21) and Eq. (22), but for the 
present case with the harmonic oscillator renormalizcd. 
These relations read 



LrY 



HMlu% 



x 



—Y ) , (62) 



Py 



L- X x = - 



eB z 



eB z 



X 



hmu;l \ Px c ^ 



(63) 



hmujy 



eB z 



Py 



Making use of the relations Eq. (62) and Eq. (63), and 
also with the effective spin-orbit interaction expressed in 



the center-of-mass and relative coordinates, see Eq. (9), 
the spin-Hamiltonian H s takes the form 



H„ = 



1 



1 



1 



1 » 12 

i,2 I v y 

v 



for the case of a perpendicular magnetic field and 

cos 2 $ sin 2 $ 



(64) 



El 



1 - 



b 2 



bl 



ctIoJ, (65) 



for the case of an in-plane magnetic field which makes 
an angle $ with the inter-dot distance direction. The y 
is along the in-plane direction perpendicular to the in- 
plane magnetic field. We see that the spin Hamiltonian 
depends on the Coulomb interaction part via the differ- 
ence between the inverse of the renormalized frequencies 
w Xi y and the bare one u>o- As expected, when there is 
no renormalizations of the bare frequencies (no Coulomb 
interaction) the interaction vanishes. Referring again to 
the effective spin Hamiltonian from Eq. (18), we 

obtain for the coupling J e g for arbitrary magnetic field 
orientations 

Jeff ' = :E " I^(a^) J ( a °/ A ' aB / A )' ( 66 ) 



where 



I(a /X,a B /X) 



?)<■ 



cos 2 6 cos 2 $ + sin 2 <f>) 



1 - 72 ) (cos 2 6 sin 2 $ + cos 2 $) . (67) 



One can see from Fig. 5b that for in-plane magnetic 
fields one obtains quite large values for I in the two lim- 
iting cases $ = and $ = tt/2. Changing the mag- 
netic field orientation in-plane one can tune the coupling 
strength J e g from negative to positive values, i.e. from 
ferromagnetic to antifcrromagnctic regime, and make it 
vanish for the angle (for in-plane magnetic field) 



$ = arcsin (b yy J(l & 2 )) . 



(68) 



In the case of a perpendicular magnetic field, cf. Fig. 5a, 
we see that the coupling shows a non-monotonic behavior 
as a function of distance ao, and, moreover, J e g vanishes 
for some given distance, which for \/cib — 5 is about 
a /A w 2.5. This could be used to tune J e g on and off 
by changing the distance between the dots. 

Next, we consider the case of very elliptic dots, with the 
bare oscillator frequencies wo^.oy corresponding to the x 
and y directions, respectively, such that uj 0x <S u>o y . The 
spin-spin coupling becomes in this limiting case 



H s = 



m* 2 X 



SO^x 



1 



cos 2 $ alal+0({u 0x /ujoy) 



(69) 
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where both perpendicular ($ = 0) and in-plane magnetic 
fields are contained. We see that the problem becomes 
effectively ID with an Ising-type spin-spin coupling, sim- 
ilar to the case of small Coulomb coupling studied in 
the previous section. We mention that for finite ratio of 
the two bare frequencies, u>o x /ujoy, the interaction can be 
varied by changing this ratio, the angle of cancellation 
defined in Eq. (68) varying as well. 

The behavior displayed in Fig. 5 can be understood 
as follows. The spin-spin coupling is directly related to 
the deformation of the charge distributions in the two 
dots as a consequence of the strong Coulomb interac- 
tion (S 3> 1). Thus, the stronger the deformation is, the 
stronger the spin-spin coupling becomes. Or, in our case, 
the stronger the deviation of the renormalized orbital fre- 
quencies Lo x>y from the bare one u>q is, the stronger the 
coupling becomes, see Eqs. (64) and (65). While the 
x component of the spin-spin coupling is bounded be- 
cause the inverse of w x tends to zero as the Coulomb 
interaction strength 5 increases, the y component of this 
coupling is unbounded since the inverse of to y can grow in- 
definitely. Consequently, the y component will dominate 
the x component for large Coulomb strength and small 
inter-dot distance do- However, the situation is reversed 
in the large distance limit, since ui x increases faster than 
Lo y decreases as seen from Eq. (55). These opposite lim- 
its lead to the non-monotonic behavior depicted in Fig. 
5a (perpendicular magnetic field). 

We mention that the large distance limit of Eq. (66) 
converges to the large distance result obtained in the pre- 
vious section, Eq. (38). However, it does not converge to 
the results of the previous section in the case of small 
distance [Eq. (40)], when going from <53>lto<5<§;l,a 
crossover description (6 ~ 1) being needed in this situ- 
ation. Phrased differently, tuning the spin-spin coupling 
J e ff from strong (S ^> 1) to small (S <C 1) Coulomb 
interaction regime by varying the inter-dot distance re- 
produces the corresponding 5 < 1 result in Eq. (38), 
while by varying the ratio A/as does not reproduce the 
corresponding S <C 1 limit, i.e. Eq. (40). 

Let us give now some estimates for the coupling J e g 
when an in-plane magnetic field is applied along, say, 
the x-direction. Assuming now GaAs quantum dots, and 
E z = 0.1 meV {B = 4 T), haj = 0.5 meV (X/a B w 5), 
A/Aso ~ 10 _1 . Using these numbers and taking for 
the geometric inter-dot distance a^/X w 2, we obtain 
J e ff w 10~ 7 eV. It is worth mentioning that the hyper- 
fine interaction between the electron and the collection 
of nuclei in a quantum dot (w 10 5 ) leads to similar en- 
ergy scales 43,44 . This shows that the spin-spin coupling 
derived here can be very relevant for the spin dynamics 
in the case of electrostatically coupled quantum dots and 
that it can also compete with other types of interactions. 
Considering now the case of InAs quantum dots in a mag- 
netic field along the x direction, with A50 ~ 2A w lOOnm 
and Ez/huJo = 0.1 and taking also ao/A w 2, a value of 
Jeff ~ 10~ 6 eV it is obtained. However, this is just a 
rough estimate since the spin-orbit coupling cannot be 



treated as a perturbation anymore and our approxima- 
tion, being pushed to the limit of its range of validity, 
starts to break down. 



IV. MEASUREMENT SCHEME 

In this section we propose a measurement scheme for 
the spin-spin interaction Jeff- Similar to the spin relax- 
ation experiments in Ref. 8, the left dot is monitored 
by a sensitive charge detector, such as a quantum point 
contact (QPC) or a single electron transistor (SET). We 
show the main steps of the scheme in Fig. 6. 

The first step is the initialization step shown in 
Fig. 6(a). At low temperatures, T <C Ez, a single- 
electron dot will relax to the ground state after a time 
larger than the spin relaxation time T\ ~ (0.1 — 100) ms. 
A faster spin relaxation can be induced by cotunneling 
with the lead, for which the dot can be placed closer 
to the Fermi surface for some time. In Fig. 6(a), the 
left dot is initialized in the lower Zeeman sublevel | f), 
whereas the right dot is empty. Next, the right-dot en- 
ergy is lowered below the Fermi energy and the dot is 
quickly filled with an electron in either upper or lower 
Zeeman sublevel. This is a sequential tunneling process 
and we denote its rate by T. In Fig. 6(b), both dots 
are deep below the Fermi surface and J e ff is the energy 
scale that governs a coherent evolution in the subspace 
{| Tl)) I IT)}- The two spin density matrix reads 

e(t) = ^ITT)mi + ^l*(*)X*(*)l, (70) 

where IvE'(t)) is the wave function that describes the oc- 
currence of the state | tl) in the initialization step. In the 
ideal case, \&(t)) evolves coherently due to the spin-spin 
interaction, Eq. (18) 

*(t)) = cos(4Jefft/H)\ tl) - ism(4J eff t/h)\ |t). (71) 

Here, we neglect the cotunneling and other spin re- 
laxation processes. In particular, the cotunneling rate 
T 2 /U± should be much smaller than the spin-spin cou- 
pling J e ff. Here, U± is the addition/extraction energy of 
the single-electron quantum dot in the Coulomb blockade 
valley. On the other hand the sequential tunneling rate T 
should be large enough, so that the spins have no time to 
evolve during the initialization and measurement steps. 
We summarize the required regime by the inequality 

^ « Jeg « r. (72) 

After a waiting time r, the probability of the left-dot 
electron to be in the upper Zeeman sublevel reads, 

P Ll (r) = \[l-cos(8 Jeff t/H)}. (73) 

Form the period of this function (to = irh/AJ e ff) one can 
extract the value of the coupling constant J e g. 
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FIG. 6: Scheme to measure the coupling constant J e g in quan- 
tum dots without tunnel coupling. In the initialization step 
(a), the left dot is at equilibrium with one electron in the 
lower Zeeman sublevel and the right dot is empty. At the 
start of the coherent evolution step (b), the right dot is filled 
with one electron in either upper or lower Zeeman level dur- 
ing a short time h/Y <C h/J e g, and the dots are deep in the 
Coulomb blockade valley, Y 2 /U± <C J e g. Further, the two 
spins evolve coherently due to the spin-spin interaction J e g 
during a fixed time r > h/J e g. In the read-out step (c), the 
left dot is brought up to the Fermi surface, so that the elec- 
tron can tunnel to the lead only if it is in the upper Zeeman 
sublevel. The latter event is recorded by a charge detector 
nearby the left dot. 



The measurement of the probability (r) can be per- 
formed in the same fashion as in Ref. 8. After the waiting 
time r, the left dot is brought up to the Fermi level and 
placed such that the electron can tunnel into the lead 
only from the upper Zeeman sublevel. This configura- 
tion is shown in Fig. 6(c). Tunneling of the electron out 
and refilling the quantum dot with an electron of the op- 
posite spin is monitored by the charge detector close to 
the left dot (not shown) . For each value of the waiting 
time r, the cycle of initialization, coherent evolution, and 
measurement has to be repeated many times in order to 
reach a good accuracy. 

Next, we remark that the hyperfine interaction with 
the lattice nuclei should not impede the measurement of 
J e ff as long as AJ e g <; A/y/N, where A is the atomic 
hyperfine coupling constant and N is the number of nu- 
clei in both quantum dots. Note that the ratio of J e g 
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FIG. 7: Residual oscillations in the averaged probability 
Pli(t) for values of A/y/N that exceed the spin-spin interac- 
tion strength 4J 6j gr. The period of oscillations is not affected 
by the hyperfine interaction and is given by to = Trh/4J e ff. 
With increasing the hyperfine strength A/y/N, the amplitude 
of oscillations decreases as oc y/N /A. As a function of the 
waiting time r, the envelope of oscillations decays as oc 1/y/r. 



to A/y/N for a constant A/a scales with the dot lat- 
eral size as oc A 4 , for strong Coulomb interaction, and 
as oc A 6 , for weak Coulomb interaction. Therefore, the 
regime 4J eJ y > A/y/N can be easily achieved by tak- 
ing a larger quantum dot. Furthermore, the hyperfine 
interaction with the nuclei has only the effect of reduc- 
ing the visibility of oscillations of Pn(t), and even for 
A/y/N » J e ff a small part of Pl[{t) shows oscillatory 
behavior with unchanged period, t = irh/AJ e ff. 

In Fig. 7, we plot the probability Pl[{t) averaged 
over the realizations of the hyperfine field. We choose 
A/y/N > 4J eJ gr to show that the measurement scheme 
is robust against the hyperfine field. The oscillations are 
well visible even when A/y/N is several times larger than 
4 J e ff . The averaged probability Pli(t) is obtained in the 
following way. For the subspace {| fj.), I it)} an d in the 
limit Ez 2> A/y/N, the coupling of spins to the hyper- 
fine field is given by Hgh z = \8h z {a{ — erf), where the 
hyperfine field Sh z has a Gaussian distribution with zero 
average and a variance a — A/y/N, which we take to be a 
measurable parameter that defines A^. 45 For more detail 
on the derivation of Hsh z we refer the reader to Ref. 46. 

For our description to be accurate, the time between 
subsequent cycles of initialization, coherent evolution, 
and measurement should be larger than the nuclear spin 
relaxation time. Considering the sum of Hsh z and 
in Eq. (18), we find that the probability Pl[{t) for a 
fixed value of Sh z is given by 47 



1 



Pli(t) = 



— cos ^ 



->.-/, ' Jl6J 2 eff + Shlj 



(l + Shl/ 



16J 



'■JJ 



(74) 



The averaged probability Pli(t) is then computed by 
integrating Eq. (74) over Sh z with the Gaussian weight 
factor P a {8h z ) = — 1 7 =e~ 5h ^ 2ry2 , where a = A/y/N. 
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Considering A/VN > AJ e g and r > irh/AJ e g, we find 
that the visibility of oscillations scales with A/y/N and 
t as follows 



v oc 




(75) 



Note that the scaling law v oc \/N /A is weaker than what 
one might expect naively from Eq. (74), after substituting 
there 5h z by its typical value A/\/N, which gives v ex 
N/A 2 . The reason for the weaker scaling law is the fact 
that 6h z is centered around zero and the denominator in 
Eq. (74) has nearly no effect. We find that the numerical 
results in Fig. 7 can be reproduced fairly accurately, if 
we approximate the argument of the cosine in Eq. (74) 
as follows 



1GJ eff 



Shi « 4 J eff 1 



Sh 2 



19 J 2 



(76) 



This approximation is justified in the regime r > To 
by the the minimal phase requirement, despite the fact 
that Sh z / 'Jeff may be large. With Eq. (76), it is easy to 
average Pli(t) and obtain an approximate expression, 
which is fairly accurate for r > t and asymptotically 
exact in the limit r 3> t . Wc thus obtain 



Pli(t) = 1\P-6 P (t)}, 



Cexp 



Sp(t) = 



2"V?) crfc (7! 

Ccos [2ttt/t + <A)(t)] 



[C 4 + (2Wto) 2 ] 1/4 



(77) 
(78) 
(79) 



where £ = 4J e ffV N /A, erfc(£) is the complementary er- 
ror function, and the running phase shift <po(t) is given 

by 



, , 1 /2ttt\ 
<Po{t) = -arctan I — J ■ 



(80) 



We note that Eq. (77) is exact in two limiting cases: £ 3> 
1 and t > To. In particular, in the limit £ — > oo, we 
recover Pli(t) in Eq. (74) for all values of r. In the 
opposite limit, £ <C 1, Eqs. (77)— (80) can be significantly 
simplified, yielding 



PLi(T) = \ 



f-\/^ cos ( 27rT / T o + tt/4) 



(81) 



V. DISCUSSIONS AND CONCLUSIONS 

In the entire derivation we assumed no tunneling be- 
tween the dots. Even in the presence of tunneling, when 
direct Coulomb repulsion, £/ 12 (which is just the classical 
interaction between two charge distributions) is larger 



than the exchange interaction, J exc (C/12 S> J exc), the 
theory presented here is expected to remain still valid. 
The reason is that this type of spin-spin coupling is a 
direct consequence of the deformation of the electronic 
charge distribution due to Coulomb repulsion between 
the two electrons. Since this is given by the sum of the 
direct Coulomb part U12, and the exchange part J eX c, 
the spin coupling J e g will be insensitive to exchange in 
the limit U12 S> Jexc- The point now is that with finite 
tunneling, even with the assumption that the coupling 
strength J e g is not modified by the exchange Coulomb 
interaction, the resulting Hciscnberg exchange coupling 
Jexc {Hexc = JexcS\ • S2) will start to compete with the 
electrostatically induced spin interaction. As a conse- 
quence, the spin coupling J e g will be washed out in the 
limit Jexc S> Jeff- However, since the Hciscnberg ex- 
change coupling decays with the inter-dot distance like 27 
Jexc ~ exp(— 2a 2 ,) while J e g ~ ag 3 , the electrostatic 
spin coupling will start to dominate at not very large 
distances. 

We recall that when having exchange, another type 
of spin coupling, induced by spin-orbit interaction 

comes into play, and which is proportional to the Heisen- 
berg coupling J e xc, i.e. Jf x ° c cx (0.1 - 0.01) J exc for typ- 
ical GaAs dots 27,30-34 . Thus, can be much larger 
than J e ff for large enough Heisenberg coupling, with a 
crossover from this exchange- type to the direct Coulomb- 
type coupling taking place at some inter-dot distance. 
This crossover, however, occurs before we get J eX c ~ Jeff 7 
since Jfxc i s typically 100 times smaller than Jexc- To 
give an estimate, we assume J exc ~ 10~ 5 eV for do = 1, 
which gives J exc ~ 3.5 • 10~ 8 eV for a = 2, implying that 

J exc J eff • 

Another important issue is the effect of screening in- 
duced by the surrounding electrons in the 2DEG and the 
metallic gates. As is well-known, the screening effect be- 
tween two charges becomes important for distances ex- 
ceeding the screening length A scr ~ XFermi (Fermi liq- 
uids). However, the screening of bare Coulomb inter- 
action depends strongly on the dimensionality. In 3D 
the effect of screening is to induce an exponential de- 
cay of the bare Coulomb interaction 48 , with the decay 
parameter A scr , while in 2D the decay follows a power 



law (<• 



in the large distance limit ) 



49-52 



with \ s 



being the relevant length scale. For GaAs, the screen- 
ing length is around A scr <~ XFermi ~ 50 nm. Moreover, 
additional screening is introduced by the electrodes to 
gate the dots, due to their metallic character. The fi- 
nite screening implies then that our theory in fact over- 
estimates the strength of the electrostatically induced 
spin coupling J e ff for distances exceeding this screening 
length and the results obtained here become just an up- 
per bound on J e g for this limit . 

Being highly controllable, the coupling J e g could be 
used to perform two qubit gates for the realization of 
quantum computers with electron spins, like proposed in 
Ref. 4. The switching times range between rather slow 
(~ 10 /is in GaAs) and reasonably fast (~ 50 ns in InAs). 
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When making use of the standard exchange coupling 4 
for switching (with typical switching times of 100 ps in 
GaAs) the electrostatically induced spin-coupling found 
here can lead to gate errors. However, this effect can 
be controlled by choosing the magnetic field direction 
or strength and/or the interdot distance such that J e g 
becomes negligibly small (see Eq. (67)). 

Finally, an important question is how orbital fluctua- 
tions (for example of the confining energy hu>o) mediated 
via spin-orbit coupling lead to fluctuations in the cou- 
pling J e ff and thus eventually to spin decoherence. The 
relation between the orbital dcphasing time (which is as- 
sumed to be known) and the decoherence induced by the 
spin coupling J e g reads 53 



SJ, 



(82) 



where is the orbital dcphasing time and the corre- 
sponding spin decoherence time. We also assumed only 
fluctuations in the dot size, which means fluctuations 
of the confinement frequency ljq. Assuming an orbital 
dcphasing time w 1 ns and also the limiting case of 
touching dots with the same GaAs parameters as be- 
fore we obtain a spin decoherence time (lower bound) 
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-3 



We can conclude then that the incoherent 



part due to this type of coupling is negligible compared 
with other types of decoherence mechanisms, e.g. induced 
by the hyperfine interaction 43 ' 44,53 . 

To conclude, we have derived an effective spin-spin in- 
teraction between two electrons localized in two quan- 



tum dots, spatially separated, induced by the direct 
Coulomb interaction and mediated by the spin-orbit cou- 
pling. This interaction was found to have the form of an 
anisotropic XY interaction and to be proportional to the 
Zeeman energy. The spin-spin coupling was studied both 
in the weak and strong Coulomb interaction limits and 
for different magnetic field orientations and strengths. 
The important features are the non-monotonic behavior 
of this spin interaction for some magnetic field orienta- 
tions, together with a vanishing of this interaction for 
particular inter-dot distances. This effect can be used to 
manipulate the spin-spin interaction in electrostatically 
coupled quantum dots by tuning the inter-dot distance. 
We proposed a measurement setup which allows one to 
access this spin-spin coupling experimentally 
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APPENDIX A: J 5ij FOR ARBITRARY B-FIELDS 

In this Appendix we give explicit formulas for the cou- 
plings Jj.y for an arbitrary magnetic field orientation 
B = B(cos$sin^,sin$sin0, cos#) and for both Rashba 
and Dresselhaus spin-orbit couplings present. These are 
obtained by diagonalizing the tensor M, which leads to 



J i-,v = g + C ^ cos2 6 + ( Cl cos2 ^ + C ' 2 sin2 ^ ~ ° 3 sin 2< ^ sin2 9 



± yj ((Ci + C 2 ) cos 2 + sin 2 0(Ci cos 2 + C 2 sin 2 <j) - C 3 sin 20)) 2 - 4{dC 2 ~ Cf) cos 2 0} , (Al) 



with (f> = $ — 7 (the angle 7 is defined after Eq. (9)). 
The functions Cj (i = 1,2,3) can be expressed in terms 
of C ai b 2 = (0|[L^ 1 ai,6 2 ]|0), a,b = x,y, i.e. 

Ci = -^- (sin 2 7 C XlX2 + cos 2 7 C yiy2 + sin 2 7 C XlV2 ) 

(A2) 



° 3 = 2 A* A ( sin2 ^( C ^2 -Cyivz) ~ cos 27 C xiy2 ) . 

(A4) 

These functions can be identified very easily from our 
formulas derived in the paper. For example, for the case 
considered in Eq. (27) (weak Coulomb coupling regime) 
we get 



C z = TT ( cos2 7 C XlX2 + sin 2 7 C yiy2 - sin 27 C XlV2 ) 



(A3) 



c 



AE y r , 



c 



XlX2 ~ m* 2 A 2 W 4 ' VlV2 ~ ™*2\2,.,4> ^*iV2 







m* 2 \ 2 u>, 



C x -. ii 0. 







(A5) 
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